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Abstract: Probablistic solutions of the so called Schrodinger boundary data 
problem provide for a unique Markovian interpolation between any two strictly 
positive probability densities designed to form the input-output statistics data 
for a certain dynamical process taking place in a finite-time interval. The key 
problem is to select the jointly continuous in all variables positive semigroup kernel, 
appropriate for the phenomenological (physical) situation. 

The issue of deriving a microscopic dynamics from the input-output statistics 
data (analyzed in terms of densities) was addressed, as the Schrodinger problem 
of a probabilistic interpolation, in a number of publications [1]- [6]. We shall con- 
sider Markovian propagation scenarios so remaining within the well established 
framework, where for any two Borel sets A, B C R on which the strictly posi- 
tive boundary densities p(x, 0) and p{x,T) are defined, the transition probability 
m(A, B) from the set A to the set B in the time interval T > has a bi-variate 
density given in a specific factorized form: m{x,y) = f(x)k(x,0,y,T)g(y), with 
marginals: 

m(x, y)dy = p(x, 0) , J m(x, y)dx = p(y, T) (1) 

Here, f[x), g(y) are the a priori unknown functions, to come out as solutions of the 
integral system of equations (1), provided that in addition to the density boundary 
data we have in hands any strictly positive, continuous in space variables function 
k(x,0,y,T). Additionally, we impose a restriction that k(x,0,y,T) represents a 
certain strongly continuous dynamical semigroup kernel, while given at the time 
interval borders: it secures the Markov property of the sought for stochastic pro- 
cess. 

It is the major mathematical discovery [2] that, without the semigroup as- 
sumption but with the prescribed, nonzero boundary data p(x, 0), p{y, T) and with 
the strictly positive continuous function k(y, 0, x, T), the system (1) of integral 
equations admits a unique solution in terms of two nonzero, locally integrable 
functions f(x),g(y) of the same sign (positive, everything is up to a multiplicative 
constant) . 

If k(y, 0, x, T) is a particular, confined to the time interval endpoints, form 
of a concrete semigroup kernel k(y, s, x, t), < s < t < T ', then there exists a 
transition density: 

P(y, s, x, t) = k(y, s, x, t) d , X ' * (2) 
%,s) 
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defined in terms of functions 



9(x,t) = dyk(x,t,y,T)g(y), 9*(y,s) = / dxk(x,0,y,s)f(x) (3) 



which implements a consistent Markovian propagation of the probability density 
p(x,t) — 6(x,t)0*(x,t) between its boundary versions, according to the standard 
transport recipe: p(x,t) = J p(y, s,x,t)p(y, s)dy. For a given semigroup which is 
characterized by its generator (Hamiltonian) , the kernel k(y, s, x, t) and the emerg- 
ing transition probability density p(y, s, x, t) are unique in view of the uniqueness 
of solutions f(x),g(y) of (1). For Markov processes, the knowledge of the tran- 
sition probability density p(y, s, x, t) for all intermediate times < s < t < T 
suffices for the derivation of all other relevant characteristics. 

In the framework of the Schrodinger problem the choice of the integral kernel 
k(y,0, x,T) is arbitrary, except for the strict positivity and continuity demand. It 
is thus rather natural to ask for the most general stochastic interpolation, that is 
admitted under the above premises. 

Clearly, the familiar strictly positive (Feynman-Kac) semigroup kernels gen- 
erated by Laplacians plus suitable potentials are very special examples in a sur- 
prisingly rich encompassing family. Indeed, the concept of the " free noise" , nor- 
mally characterized by a Gaussian probability distribution appropriate to a Wiener 
process, can be extended to all infinitely divisible probability laws via the Levy- 
Khintchine formula [7]. It expands our framework from continuous diffusion pro- 
cesses to jump or combined diffusion-jump propagation scenarios. All such (Levy) 
processes are associated with strictly positive dynamical semigroup kernels, and all 
of them give rise to Markov solutions of the Schrodinger stochastic interpolation 
problem (l)-(3). 

At this point, let us remark that apart from the wealth of physical phenomena 
described in terms of Gaussian stochastic processes, there is a number of physical 
problems where the Gaussian tool-box proves to be insufficient to provide satisfac- 
tory probabilistic explanations. Non-Gaussian Levy processes naturally appear 
in the study of transient random walks when long-tailed distributions arise [8,9]. 
They are also found necessary to analyze fractal random walks [10], intermittency 
phenomena, anomalous diffusions, and turbulence at high Reynolds numbers [8,11]. 
On the other hand, our formulation of the Schrodinger interpolating dynamics can 
be regarded as a straightforward inversion of the well developed programme of 
studying dynamical systems (chaotic included) in terms of densities, [12,13]. 

Let us consider Hamiltonians (semigroup generators) of the form H = F(p), 
where p = — iV stands for the momentum operator and for — oo < k < +oo, 
F = F(k) is a real valued, bounded from below, locally integrable function. Then, 
exp(-tH) = exp[—tF(k)]dE(k), t > 0, where dE(k) is the spectral measure 
of p. We simplify further discussion by considering processes in one spatial di- 
mension. Because (E{k)f)(x) = -j== exp(ipx)f(p)dp, where / is the Fourier 




transform of /, we learn that 
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exp(-tF(k))exp(ikx)f(k)dk = [exp(-tF(p))f(p)] v (x) (4) 



where the superscript V denotes the inverse Fourier transform. 

Let us set k t — —!=={exp(—tF(p)] v , then the action of exp(-tH) can be given 
in terms of a convolution: exp{—tH)f — f * k t , where (/ * g)(x) :— f R g{x — 
z)f(z)dz. 

We are interested in those F(p) which give rise to positivity preserving semi- 
groups: if F(p) satisfies the celebrated Levy-Khintchine formula, then k t is a 
positive measure for all t > 0. The most general case refers to a contribution 
from three types of processes: deterministic, Gaussian, and an exclusively jump 
process. Let us concentrate on the integral part of the Levy-Khintchine formula, 
which is responsible for arbitrary stochastic jump features: 

F(p) = [ + °°[exp(ipy) 1 - -^Hdy) (5) 

where v{dy) stands for the so-called Levy measure. The disregarded Gaussian 
contribution would read F(p) — p 2 /2; cf. Refs. [1-6] for an exhaustive discussion 
of related topics based on the concept of the conditional Wiener measure. 

There are not many explicit examples (analytic formulas for probability den- 
sities involved) for processes governed by (5), except possibly for the so called 
stable probability laws. The best known example is the classic Cauchy density. 
Let us focus our attention on two selected choices for the characteristic expo- 
nent F(p), namely: Fq(p) = \p\ which is the Cauchy process generator, and 
Fm(p) = \lv 2 +m 2 — m,m > 0. Here, we have chosen suitable units so as to 
eliminate inessential parameters. The latter exponent is ano ther form of the fa- 
miliar classical relativistic Hamiltonian, better known as \Jm 2 c 4 + c 2 p 2 — mc 2 
where c is the velocity of light. The respective semigroup generators Ho, H m are 
pseudodifferential operators. The associated kernels fc™ in view of the "free 
noise" restriction (no potentials at the moment) arc transition densities of the jump 
(Levy) processes with intensities regulated by the corresponding Levy measures 
^o(rfy), v m {dy)- The affiliated Markov processes solving the Schrodinger problem 
(l)-(3) immediately follow. It is instructive to notice that like in case of more tra- 
ditional Gaussian derivations [4], the identities 9(x, t) = 1, 9*(x, i) := p(x, t) imply 
the pseudodifferential analog of the Fokker-Planck equation. It is a consequence of 
[exp(— tH)~p\(x) = p(x, t) and of the identification F(p — > — iV) := H . For example 
there holds: F (p) =^> dip(x, i) = — |V|p(x, t). This evolution rule gives rise to the 
Cauchy process with its long-tailed probability density rho(x, t) — ^ ^i\ x i and the 
transition probability density (e.g. the semigroup kernel in this free propagation 
case) of the same functional form with x — > x — y and t — > t — s. Let us emphasize 
that the existence and uniqueness of solutions proof for the Schrodinger problem 
extends to all cases governed by the infinitely divisible probability laws, and can 
be generalized to encompass the additive perturbations by physical potentials (in 
analogy with the familiar Feynman-Kac formula). 
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Our semigroups are holomorphic, hence we can replace the time parameter t 
by a complex one a = t + is, t > so that exp(-aH) — J R exp(—aF(k))dE(k). 
Its action is defined by [exp(—aH)]f = [( fexp{— <tF)] v = f * k a . Here, the ker- 
nel reads k a = -j=[exp(— <?F)] V . Since H is sclfadjoint, the limit t j leaves 
us with the unitary group exp(—isH), acting in the same way: [exp{— isH)]f = 
[fexp(— isF)] v , except that now k is := -J=[exp(— isF)] v in general is not a mea- 
sure. In view of unitarity, the unit ball in L 2 is an invariant of the dynamics. Hence 
density measures can be associated with solutions of the Schrodinger pseudodifcr- 
ential equations: F (p) id t ip(x,t) = \ V\4>(x, t) or F m (p) id t ip(x,t) = [y/— A + m 2 — m]ip(x,i), 
if provided with the appropriate initial data functions ip(x,0). Let us point 
out that we know in detail how the analytic continuation in time of the Lapla- 
cian generated holomorpic semigroup induces a mapping to diffusion processes 
of a quantum mechanical proveninence (since the standard Schrodinger equation 
id t ip(x,t) — —Ai[)(x,t) is involved). 

All that ultimately submits the unitary (quantum) Schrodinger picture dy- 
namics, with quite a variety of admissible semigroup generators to be used instead 
of the traditional Laplacian, to the stochastic analysis in the framework of the 
Schrodinger (again) boundary data problem. The natural question to be answered 
is: what are the stochastic processes consistent with the probability measure dy- 
namics p(x,t) = \ip(x,t)\ 2 determined by pseudodiffcrential equations, eventually 
in the presence of external force fields ? The answer to this and related ques- 
tions of the more pedestrian, nonequlibrium statistical physics origin can be found 
elswhere, [6,13,14]. 
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